clear all
clc
[x,f_min]=fminbnd('-(4-2*x)^2*x',0,2)
